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In this paper we investigate the emergence of time-periodic and and time-quasiperiodic (some- 
times infinitely long lived and sometimes very long lived or metastable) solutions of discrete nonlinear 
wave equations: discrete sine Gordon, discrete 0* and discrete nonlinear Schrodinger. The solutions 
we consider are periodic oscillations on a kink or standing wave breather background. The origin of 
these oscillations is the presence of internal modes, associated with the static ground state. Some of 
these modes are associated with the breaking of translational invariance, in going from a spatially 
continuous to a spatially discrete system. Others are associated with discrete modes which bifurcate 
from the continuous spectrum. It is also possible that such modes exist in the continuum limit and 
persist in the discrete case. The regimes of existence, stability and metastability of states as the 
lattice spacing is varied are investigated analytically and numerically. A consequence of our analysis 
is a class of spatially localized, time quasiperiodic solutions of the discrete nonlinear Schrodinger 
equation. We demonstrate, however, that this class of quasiperiodic solution is rather special and 
that its natural generalizations yield only metastable quasiperiodic solutions. 
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I. INTRODUCTION 

In the past two decades, there has been an explosion 
of interest in the study of solitary waves in spatially dis- 
crete nonlinear systems. Prominent examples, where one 
expects important phenomena to be intrinsically associ- 
ated with spatial discreteness, are: the motion of dis- 
locations in solid-state physics the transmission of 
kink-like bits of information in Josephson junction arrays 
1^ , the propagation of pulses in nonlinear optical waveg- 
uides |^,|^ (where they present great potential for appli- 
cations such as all-optical switching, steering and demul- 
tiplexing see e.g., Q| and references therein), or even in 
chemical and biological jsj (and references therein) 
applications. Recent experimental work confirming the 
relevance of such models has appeared in @J|,|l^ . 

The models arising in many of these studies, governing 
the behavior of a field Un{t) depending on a continuous 
time variable t and discrete spatial variable n, are: 
the discrete SG equation (DSG) 



Un = Un+l + Un-l - 2u„ - — sinM„, 



(1) 



the discrete (j)"^ equation {Dcj)'^) 



Un = Un+l + Un-1 - 2u„ + ^ ("« - U^), (2) 

and the discrete nonlinear Schrodinger equation (DNLS) 

iiln = -C(m„+1 + Un-l - 2Un) - |u„|^U„. (3) 

We have used the standard notation where d = 1/Ax is 
the inverse spatial lattice spacing and C — d^ — (Aa;)~^ 
is the coupling constant. 

An important phenomenon associated with spatial dis- 
creteness is the role of the spatial lattice sites as "de- 
fects" to which propagating states can be "trapped" or 
"pinned" ||ll|-|l^. At a structural level, this is a result 
of the breaking of translational invariance. Indeed, the 
spatially continuous variants are Lorentz invariant (SG, 
(f>*) and Galilean invariant (NLS); a stationary coherent 
structure can be "boosted", i.e., transformed, into one 
traveling with a constant velocity. 

The main coherent structures of interest in these equa- 
tions (sec i.e., ||l3|- |l5| , |l9| -^ ) are kink- like excitations in 
DSG and Dcf)* and time periodic, exponentially localized 
in space pulse like excitations in DNLS (solitary waves). 
The latter are frequently referred to as discrete breathers. 

In the analytical and numerical work on these discrete 
systems Il^,0,0,|l3 two main steady states have been 
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identified; there are actually many other steady states 
1^,^. One is centered on a lattice site and one is 
centered between consecutive lattice sites. In the case 
of DSG and D(j>'^ kinks, the stable configurations (en- 
ergy minimizers) are those centered between sites, while 
for the breathers of DNLS, stable configurations (con- 
strained energy minimizers; the constraints are imposed 
by the additional conserved quantities such as the P 
norm of the solution ) are those pulses which are cen- 
tered at a lattice site. 

Until recently, it was believed that an initial condi- 
tion, which for the translation invariant systems leads to 
a coherent structure propagating to infinity at a uniform 
speed, would in these discrete systems result in a state 
radiating energy which would eventually get trapped by 
the "Peirels -Nabarro" potential associated with the 
lattice and subsequently approach a time independent 
state. However, it was established in Q that, depending 
on the parameter regime (lattice spacing) , the asymptotic 
state may, in fact, be time-periodic. It may also consist 
of a time quasiperiodic state which slowly decays to a 
time-periodic state. The rate at which the system re- 
laxes to its asymptotic state is controlled by resonances 
of "localized internal modes" with radiation modes. In- 
ternal modes are spatially localized modes associated 
with the linearized operator about the coherent struc- 
ture with corresponding discrete eigenfrequencies (point 
spectrum) , which are not related to an underlying group 
invariance (zero modes). They give rise to time periodic 
solutions of the linearized evolution equation, of the form 
(j)':^ exp{iujt) . Radiation modes are those associated with 
the continuous spectrum of the linearization. 

In a recent pa per [l^ , using a methodology developed 
in jl^] (see also [|l7|l ) , we gave a systematic approach for 
deriving details of the decay to an asymptotic state. We 
now briefly describe these results. Suppose oji and uj2 are 
two internal mode frequencies associated with a coherent 
structure of (|l]), (||) or (||). Consider initial conditions 
which are a small perturbation of the coherent structure 
and that the perturbation (of 0{e)) consists of compo- 
nents in the direction of the internal modes, (f)'^^ and cf)'^^ . 
As time evolves, nonlinearity generates various linear 
combinations of these frequencies so we investigate where 
such combinations land relative to the so-called phonon 
band (PB). Suppose that a combination of frequencies 
kuji + luj2 falls inside the PB with -f |;| = n. For DSG 
and D(f)*, PB is, respectively, the purely imaginary in- 
terval ±i[l/d,^/ITl/d?],±i[V2/d,^/t+2/¥] and for 
DNLS it is ±i[A, A -I- 4C], where A is the time- frequency 
of the ground state discrete solitary wave. In this case 
we showed in that the time evolution of the inter- 
nal mode powers (squared amplitudes), P = P^-^ and 
Q = Qi^2 is controlled by a system whose normal form is 



Pt^- e^(\^\+\^\-^^T,.MP^Q' (4) 

Qt = - £'^l'l+l'l-'^r2;fe,P^Q', (5) 

where /'"'^^ — {k, I : koji + Ilu2 G PB} is the set of 
resonant frequency combinations. The "damping coeffi- 
cients" Ti-ki, ^2;ki can be computed [T^ ], are always non- 
negative and are generically strictly positive. From (^- 
(^, one concludes that the decay rate is controlled by the 
appropriate harmonic e.g., if 2wi G PB, then the (long 
time asymptotic) decay of P goes as t^^^^, if Swi S PB, 
it goes as i^^^^ and so on. As the parameters of the 
problem (such as the lattice spacing) vary, the set of res- 
onances may vary and therefore so will the normal 
form equations and the nature of the asymptotic state 
dynamics; sec [ p^ for examples of such transitions. 

II. ASYMPTOTIC STATES: STATIONARY AND 
PERIODIC, STABLE AND METASTABLE 

What are the possible asymptotic states ? Consider 
the discrete nonlinear wave equations (|l]) and (||) . There 
is, of course, the ground state kink but one may ask 
whether there are additional time-periodic or even quasi- 
periodic states associated with the time-periodic internal 
modes., e.g., solutions of the nonlinear equation of the 
form: Kgs + I]t^eni„te™ai exp(za;i)-h C.C. plus higher 
order corrections. Here, Kgs denotes the ground state 
or least energy kink and ^internal denotes the set of lo- 
calized internal mode frequencies. In [ p^ we prove the 
existence of time-periodic solutions (only one contribu- 
tion to the above sum is present) in certain parameter 
regimes. No quasi-periodic solutions exist since for any 
internal mode frequencies loi^lo2^ there always exist inte- 
gers fci 7^ 0, fc2 7^ such that kiijJi^k2'jJ2 G PB. There are 
however metastable, very long-lived states of this type. 

For DSG, depending on parameter regime, we have 
between one and three possible asymptotic states from 
among Kgs (the ground state kink), gW and eW , each 
with their own basin of attraction. gW and eW are 
time-periodic solutions associated with perturbed trans- 
lational modes (i.e., bifurcating from the origin of the 
spectral plane) and edge (bifurcating from the edge of 
the continuous spectrum) internal modes (Theorem 5.1 
of |l|). See Table I. 

For Dcf)^ , we have between one and three from among 
Kgs, gW, W and eW. The additional periodic solution 
W , is associated with the additional shape mode present 
in the point spectrum of the cf)'^ model. See Table II. 

The above periodic solutions are localized time- 
periodic excitations on a kink (spatially inhomogeneous) 
background, "breathers on kinks" (BOK's); see Fig. 1. 
Their dynamic stability in various regimes of the discrete- 
ness parameter d was numerically investigated in [p3[. 
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We now turn to DNLS. Seek solutions of the form 

Unit) = exp(iAt)V'„. (6) 

This gives rise to a system of nonlinear algebraic equa- 
tions for {"ijjn} 

- Aljjri = -C(V'n+l + 1pn-l " S^-n) - IV'nPV'n- (7) 

A variational approach to the questions of existence, exci- 
tation thresholds and nonlinear Lyapunov stability of the 
ground states of is presented in jl^. Solutions can 
also be constructed by iterative methods implemented 
numerically in . Linear stability analysis proceeds 

by seeking a solution of the form 



Unit) = exp(iAt)(-0„ + Vnit)) 



(8) 



where w„ = a„ exp(— icijt) + 6„exp(ia;t) denotes the per- 
turbation about the ground state. This leads to the spec- 
tral problem 

ujan = -CA2a„ - 2|V'npa„ + Aa„ - ijjn'^b^, (9) 
-ubn = -CA2bn ~ 2|V'«P6„ + A6„ - iZ-n^a* . (10) 

The analysis of the modes for different values of the rel- 
evant parameter A/C is given in [|l9|. There it is found 
for the stable (site-centered) pulses of DNLS that there 
are two modes, a perturbed translational internal mode 
(the "pinning" mode) and an edge mode (the "breath- 
ing" mode). The picture for DNLS is derived in 
and is analogous to that described above and in for 
equations (^-(||). In fig. 2, the relevant (localized) eigen- 
mode frequencies are given and compared with those of 
the PB. The analysis of shows that the methodology 
of ijl^ can be applied in this case, even though the basis 
of eigenfunctions is not orthonormal (the breathing mode 
has non-vanishing overlap with the rotational symmetry 
eigenmodes with a; = 0). However, a pseudo-scalar prod- 
uct can be used jl9| in the spirit of to perform the 
same analysis. Johansson and Aubry do this in two al- 
ternative ways (higher order perturbation theory and a 
conservation law approach) but nonetheless reach a spe- 
cial case of the result of proving that when the p-th 
harmonic of a mode is inside PB, then the amplitude of 
the mode decays as |a(<)| ^ ^-1/(2^-2) asymptotically, as 
is suggested also by eqns. (^)-(|^). However, in the DNLS 
case, they found, additionally, that the existence of the 
additional two-dimensional manifold of eigenmodes with 
zero frequency, in order to avoid unphysical divergences 
of higher order perturbations (due to projections on this 
manifold), causes an increase of the breather mode fre- 
quency according to At ~ |apP. 

The methods of jl^ can be used to extend the picture 
developed in jiot . Using the information, displayed in 
figure 2, on the location of linear combinations of inter- 
nal mode frequencies associated with the ground state of 
equation (^, relative to the phonon band PB, we can 
conclude: 



• For any value of A/C for which the edge mode ex- 
ists (0 < A/C < 1.7), initial conditions close to it 
decay as t^^/^ since its second harmonic is always 
inside PB. 



• The perturbed translational mode (which disap- 
pears around A/C ^ 1.1 as is noted in ||l^), for 
0.48 < A/C < 1.1, 2wg is inside PB, hence the 
decay in this range goes as t~^/^. However, for 
0.3997 < A/C < 0.48, the third harmonic resonates 
with PB and hence the decay goes as t~^l'^ . Subse- 
quently for A/C e [0.3576,0.3997] and t-^l^ 
for A/C e [0.3304,0.3576] are implied by principle 
resonance occurring for higher harmonics. 



Therefore, for any value of the coupling there are har- 
monics of the individual modes that are in the phonon 
band and hence cause decay of the internal mode os- 
cillations. Hence, no genuinely periodic solutions along 
eigendirections can survive in this case (in the frame co- 
rotating with the discrete breather frequency A). The 
decay however can be extremely slow as we now explain. 
The perturbed translational mode, the neutral oscillatory 
mode to which the translation zero mode is perturbed for 
Ax > 0, approaches zero very fast for small A/C ^ Ax, 
the continuum limit. This approach has been shown to 
be exponentially fast as Aa; J, (w^ ~ exp(— tt^/Ax)) in 
p^ , p5[ . Therefore, as seen in figure 2 and can be derived 
from the above estimates, for very small values of A/C, 
only very high harmonics will fall inside PB. By the anal- 
ysis of this corresponds to very slow time decay rates 
of the internal mode oscillations, which may be hard to 
detect numerically. An example of this sort is given in fig- 
ure 3. The numerical experiment is performed for DNLS, 
with K/C — 0.370 and for the perturbed translational 
mode of w/C = 0.102, the 4-th harmonic is the leading 
one in the phonon band. The initial condition is taken 
to be the exact breather profile with a small perturba- 
tion along the perturbed translational mode eigendirec- 
tion. The background oscillations are clearly discernible 
as time-dependent side wobbles in the evolving profile. 
Our simulations show very slow decay of these oscilla- 
tions, in accordance with the theoretical prediction. The 
analysis of ||l^ implies that these modes will decay to the 
ground state standing wave breather but only over very 
long time scales. Modes such as the one that is observed 
in figure 3 are very long-lived time-periodic and spatially 
localized. In view of (^) they correspond to slowly de- 
caying (metastable) quasi-periodic oscillations of DNLS 
(which are periodic in the rotating frame). The periodic 
oscillation is mounted on a breather background. Hence 
these modes are (metastable) Breathers On Breathers 
(BOB's). We note that there are other constructions of 
quasiperiodic modes in discrete systems. Such attempts 
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(see, for example, [^) exploited continuation from the 
anti-continuum limit of two oscillators with different fre- 
quencies. Also, for polaron type models such as the Hol- 
stein model, a proof of existence of quasiperiodic solu- 
tions was give in . Our analysis above is not restricted 
to the regime where the sites are weakly coupled. How- 
ever, it should be noted that very close to the continuum 
limit, logarithmic effects similar to the ones described in 
for generalized continuum NLS equations may be- 
come relevant. The work presented herein mainly fo- 
cuses in the regimes of large or moderate discreteness 
and the mechanisms described will still be present even 
for very small discreteness. The work of however, 
raises the interesting question as to whether there is a 
crossover point beyond which (in approaching the con- 
tinuum limit) logarithmic effects become dominant over 
power law ones. Such issues will be addressed in a future 
publication. 



• More interestingly, as figure 4 suggests, the gaps be- 
tween the above coupling parameter intervals con- 
stitute parameter regimes for which no multiples 
of u>i, lie inside PB. Hence an initial condition con- 
sisting of the exact mode plus a perturbation along 
the LUi, eigendirection encounters no evident decay 
mechanism. Indeed, the proof of Theorem 5.1 of 
[ p^ can be carried through in this case to obtain 
the existence of time-periodic solutions of the non- 
linear perturbation equation governing Vn{t). In 
this way, we have constructed genuine time quasi- 
periodic solutions (BOB's) of DNLS of the form (|), 
where Vn{t) is time periodic. This is demonstrated 
in figure 5, where we display the results of a nu- 
merical experiment for C = 0.02, A = 1, where the 
internal mode frequency is = 0.294. We observe 
the very regular oscillations on top of the breather 
amplitude due to the BOB's quasi-periodicity. 



III. INFINITELY LONG-LIVED QUASIPERIODIC 
SOLUTIONS 

Are there time quasi-periodic states which are non- 
decaying? We now display a family of (non-decaying) 
time quasi-periodic and spatially localized BOB states. 
Our point of departure will be the modes introduced 
in [||, which we call twisted localized modes (TLM's). 
These modes are solutions of DNLS of the form (||) with 
{4'n} = (■ ■ • , a, 1, —1, —a, . . . ) with < a << 1. They 
can be viewed as soliton-antisoliton bound pairs which 
can exist due to discreteness There is no con- 

tinuum analog since in the continuum case a soliton- 
antisoliton initial condition naturally repels. 

Henceforth, for simplicity we will fix A = 1 and vary 
C only. We construct these modes numerically by solv- 
ing equation (^ and subsequently solve numerically eqs. 
(Pl)-(p^) to determine their linear stability. We find that 

• TLM's are stable for 1/C < 0.146. Their spectrum 
for this parameter range consists of two modes at 
u; = and an additional internal mode with fre- 
quency LOi, in the gap between and edge of PB. 

• For values 1/C > 0.146 a series of instabilities arise 
giving quite interesting phenomena that have been 
analyzed elsewhere |2^. 

For C in the range of stability, the location of multiples 
of the single non-zero frequency (oj^,) relative to PB is 
displayed in figure 4. It can be seen that 

• For 1/C e [59.712,67.681], 4tj,„ e PB and hence 
implying decay of the internal mode oscillation 
about the TLM with a rate: t'^^'^. For 1/C e 
[31.943,39.604], Slu^ e PB implying f-^/^ ^^^^^ 
and for 1/C G [11.989, 19.288], 2w„ e PB implying 
a t"-"^/^ decay rate. 



IV. GENERALIZATION TO MULTIPULSES 



We saw in the previous section that TLM's can be 
thought of as a pulse and an antipulse whose "backs" 
have merged. The question that then stems from this 
analogy is whether genuinely quasiperiodic solutions of 
DNLS arise naturally in such multipulse setups where 
two or more pulses have been concatenated. 

A general analysis of such configurations for discrete 
setups has been given in [ p8[ | for two pulses and general- 
ized in [3C| ] for multiple pulses. In these works it has been 
found that in order for multipulses to be stable, the near- 
est neighbor ones have to possess opposite parity, i.e., a 
phase difference of tt. Hence, we will consider only such 
configurations hereafter. 

In order to understand the possibilities for quasiperi- 
odic solutions of multipulse type, let us return to the lin- 
earization picture for a single pulse. As is well known, for 
the single pulse, there are two pairs of eigenvalues near 
the origin of the spectral plane, the rotational modes (of 
zero frequency) and the translational modes of frequency 
dependence (on h) as exp(— 7r^//i). For a concatenation 
of N pulses, there will be 2N pairs near the origin. N 
of them will be translational (stemming from the linear 
combinations of the translational cigcnmodcs of the indi- 
vidual pulses), 1 pair will be the rotational one of a; = 
and iV — 1 pairs will pertain to the so-called interaction 
eigenmodes p^j30| ], whose frequencies will, in general, 
behave as exp(— L), where L is the (typical) pulse sepa- 
ration. It should be noted that for such multipulse con- 
figurations to occur, it necessary for the pulses to have a 
separation larger than a minimal one L > Lmin- 
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For values of h < 0(1), as was found in |28|j3C| ], it is 
generically true that for L > Lmim the interaction eigen- 
values have frequencies typically lower than the trans- 
lational ones. Hence, since as we saw in section II, 
quasiperiodic solutions cannot be formed for such regimes 
of the lattice spacing, as the translational eigenvalues res- 
onate with the continuous spectrum, hence such solutions 
will also not be but metastable for the multipulsc config- 
urations as well. 

From the above argument, it can be seen that for mul- 
tipulse configurations, as was the case with the TLM's of 
section III, only strong discreteness contexts h » I pro- 
vide the possibility for quasiperiodic solutions. In such 
regimes of strong discreteness, the translational eigen- 
modes have already merged with the continuous spec- 
trum |l9|] and only the interaction eigcnmodes (and their 
harmonics) relative position to the continuous spectrum 
has to be accounted for. 

In order to understand the behavior of multipulses, a 
natural continuation parameter that we will use in this 
context is the separation between the pulses L. In the 
case of the TLM's, this separation will be considered h, 
i.e., the distance between the pulse maximum and the an- 
tipulse minimum. In this case, the form of the solution is 
shown in the top left panel of Fig. 6 and the form of the 
evolution of w/C as a function of A/C is given in the top 
right panel. In the middle two and lower two panels of 
fig. 6, similar results are shown for L = 2h, and L = 3h 
respectively. What can be clearly observed is that in 
these cases, the dependence of the internal (interaction) 
eigenmode on the coupling constant is such that for any 
value of h, there are harmonics of the internal mode fre- 
quency inside the PB. Hence, the TLM's are very special 
in the sense that the relation of their centers positions 
is such that it permits for a set of intervals of coupling 
constants to have no resonances of internal modes with 
the PB and, thus, quasiperiodic solutions. 

To show the generality of the above conclusion, we per- 
formed similar numerical experiments for 3-pulse config- 
urations (in which the middle one was of opposite phase 
than the two outer ones). The case of L = ft is once 
again shown in the top left panel of Fig. 7 and the corre- 
sponding resonance diagram of lo/C as a function of A/C 
appears in the top right panel of the same figure. The 
same plots for L — 2h and L = 3h are shown in the mid- 
dle and lower rows of Fig. 7 respectively. In this case, as 
per the general discussion presented above, there are two 
pairs of interaction eigenmodes. In the case of i = ft, it 
can be observed that their frequency increases as a func- 
tion of A/C and hence gaps appear in the presence of 
the modes' harmonics in the PB. Such gaps for instance 
arc found for ft G (4.677,5.292), ft e (6.537,7.482) for 
the first frequency (in solid line) or for ft S (4.321, 6.468) 
for the second internal mode (in dashed line) . Excitation 
of one of these internal modes in the respective inter- 
vals will give rise to genuinely quasiperiodic solutions. 



On the other hand, it can be observed that already for 
L = 2ft, once again the dependence of w/C on A/C is 
monotinically decreasing and that harmonics will always 
be present inside. The same picture of course in main- 
tained for L = 3ft and for larger interpulse distances. 

These results lead us to the conclusion that generi- 
cally quasiperiodic solutions of the TLM variety (and its 
multipulse generalization) will be present in DNLS type 
systems. However, this is a rather "delicate" result of the 
interplay between the lattice spacing ft and the interpulse 
separation L which has been found to be quite special to 
such. More general multipulse solutions of larger inter- 
pulse distance are invariably found to possess internal 
modes whose harmonics are resonant with the PB and 
hence in the latter case only metastable quasiperiodic 
configurations will be present. 



V. CONCLUSIONS 

In this paper, we have given theoretical justification as 
well as numerical evidence of the existence and asymp- 
totic stability of breather- like modes which consist of (de- 
pending on parameter regime) infinitely long lived or very 
long lived spatially localized and time-periodic oscilla- 
tions on a kink background (BOK's) or on a breather 
background (BOB's). These modes are different from 
previously reported ones in that they reside on a non- 
uniform (coherent structure) background rather than on 
a uniform background. Because of their relevance to 
asymptotic states for a wide range of initial conditions 
(possibly in experimental situations of interest) further 
study and understanding of these modes as well as an- 
alytical investigation of their stability are challenging 
questions worth pursuing. The existence of genuinely 
quasi-periodic solutions is apparently rare for nonlinear 
conservative wave equations on infinite spatial domains. 
This is in contrast to the case of nonlinear wave equa- 
tions on finite domains. The discrete case corresponds 
to finite dimensional KAM theory. The continuum case 
(that of vibrations, strings and membranes) is studied 
in ||3l|-p3[, where such solutions have been constructed. 
The non-generic nature of such quasi-periodic solutions 
was however demonstrated for discrete systems in the in- 
finite domain, and their generalization to multipulse type 
configurations was also considered. 
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FIG. 1. BOK background: Discrete SG simulation result for time evolution of an initial condition consisting of a kink plus a 
projection to the edge mode. Shown is the time evolution of the initial condition when substracting Kgs from the profile. The 
persistent oscillations of non-decreasing amplitude indicate that the BOK is indeed the asymptotic state. 
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FIG. 2. Internal Eigemodes, harmonics and combinations for the DNLS case: The translational mode and its first few 
harmonics are shown in dashed line, the edge mode and its second harmonic shown in solid line, the sum of the two frequencies 
is given by the dotted line. The band edges of the PB are shown by dash-dot lines. 
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FIG. 3. Example of a motastablc BOB: For the DNLS and A/C = 0.370 the time evolution of a BOB is shown. The top 
left panel shows the initial condition consisting of the exact breather plus a perturbation along the direction of the eigenvector 
of the translational eigenmode. The top right panel shows the time evolution of the amplitude squared of the adjacent to the 
central (no = 100 in this case) lattice sites clearly indicating the persistence of oscillations on top of the original breather. The 
bottom left panel shows the contour plot of the spatio-temporal evolution of the full BOB while the bottom right shows the 
spatio-temporal evolution of the (amplitude squared) of the BOB minus the exact breather. This also clearly indicates the 
persistence of the metastable BOB (notice that very slow decay is observed for the simulated times). 
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FIG. 4. Relative position of the internal mode frequency of the DKL solution and its p — th harmonics (solid lines) with 
respect to PB (dash-dotted lines) for different values of the coupling constant C. A = 1. 
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FIG. 5. Proof of persistence of DKL BOB mode for C — 0.02. The additional frequency of the oscillations clearly observed 
at the top is checked to agree with the predictions of linear stability theory for lj*. Shown is the spatio-temporal evolution of 
the amplitude squared (|w„(f)p) of the field (see eqn. (^). 
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FIG. 6. The generalization of the twisted modes and their respective resonance pictures. The top left panel shows the spatial 
profile of the TLM and the top right, the frequency its internal mode and its first few harmonics (given by solid line) with 
the PB (the edges of which are given by dash-dotted lines. The middle row shows the extension of the TLM for distance 
between centers equal io L — 2h (see also text); once again the left panel shows the spatial profile, while the right shows the 
internal frequency and its harmonics positions with respect to PB. The same plots are shown in the bottom row for a TLM 
with L = 3h. It is clear that in both the middle and bottom right panels the decreasing dependence of the interaction mode 
frequency (divided by C) as a function of the inverse coupling constant generically allows for the occurrence of resonances and 
hence BOB type solutions can only be metastable. A = 0.5. 
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FIG. 7. Similax plots as the ones of Fig. 6 are given here for the generahzation to three pulses. The mode with L = ft is 
given in the top left panel while the resonance picture of its two internal (interaction) eigenmodes (given by solid and dashed 
line respectively) and their harmonics with respect to to the PB (lying between the two dash-dotted lines) is given in the top 
right panel. Similar plots for the L — 2h three and the L = 3h three pulse solutions are given in the middle and bottom row. 
Notice once again that only for the L = h case there is an increasing dependence of oj/C as a function of A/C for the internal 
modes and hence there is a possibility for genuinely quasiperiodic solutions. As discussed also in the text for L > 2h, the 
decreasing nature of the above dependence generates resonances with the PB and hence becomes prohibitive for the presence 
of (other than metastable) such solutions. 
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TABLE I. SG normal form, (xjuilibria and anticipated colicrcnt structnrcs 



Regime of d 


EquilibricL 


Colierent structures 


I:d< de,de ~ 0.515 


P > 


Kgs, gW 


II: de<d< 0.565 


{(P,0),(0,Q):P>0, Q>0} 


K,s, W, gW 


III: 0.565 < d < d, 
d, ~ 0.86 


{(0,Q) :Q>0} 


Kgs, W 


IV: d, < d 


{(0,0)} 


Kgs 



TABLE II. 0* normal form, equilibria and anticipated coherent structures 



Regime of d 


Equilibria 


Coherent structures 


I: d < 0.5398 


{(0,Q) 


P > 0} 

Q > 0} 


Kgs, W,gW 


II-III: 0.5398 < d < 0.6364 


{(0,Q) 


Q > 0} 


Kgs, W 


IV: 0.6364 < d < 0.6679 


{(0,Q) 


Q > 0} 
P > 0} 


Kgs, W, gW 


V: 0.6679 < d < de 
de ~ 0.82 


{(0,Q) 


Q > 0} 


Kgs, W 


VI: rf, < d < 0.9229 


{(O.Q.O) :Q>0} 
{(().(). i?) : R > 0} 


A',,.. W. cW 


Vll: 0.9229 < d < 1.2234 
d, ~ 1.22,34 


{(0,0, ii) : P> 0} 


Kgs, eW 


Vlll: (/, < (/ 


{(().().(>)} 


A,- 
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